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Abstract 

We present a simple transformation of the formulation of the log-periodic power law formula 
of the Johansen-Ledoit-Sornette model of financial bubbles that reduces it to a function of only 
three nonlinear parameters. The transformation significantly decreases the complexity of the fit- 
ting procedure and improves its stability tremendously because the modified cost function is now 
characterized by good smooth properties with in general a single minimum in the case where the 
model is appropriate to the empirical data. We complement the approach with an additional 
subordination procedure that slaves two of the nonlinear parameters to what can be considered 
to be the most crucial nonlinear parameter, the critical time t c defined as the end of the bubble 
and the most probably time for a crash to occur. This further decreases the complexity of the 
search and provides an intuitive representation of the results of the calibration. With our proposed 
methodology, metaheuristic searches are not longer necessary and one can resort solely to rigorous 
controlled local search algorithms, leading to dramatic increase in efficiency. Empirical tests on 
the Shanghai Composite index (SSE) from January 2007 to March 2008 illustrate our findings. 

Keywords: JLS model, financial bubbles, crashes, log-periodic power law, fit method, opti- 
mization. 
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I. INTRODUCTION 



Financial crises are crippling national economies, as evidenced by the episode of the "great 
depression" that followed the great crash of October 1929 and by the "great recession" that 
followed the financial debacle that started in 2007 with the cascade of defaults of financial 
debt securities. In these two cases [l, 2|, as well as in many others [jij], one can observe 
that the crisis followed the burst of one or several bubbles, defined qualitatively as an 
exaggerated leverage in some industry sector. This observation has led policy officials to 
call for the development of methodologies aiming at diagnosing the formation of bubbles as 
early as possible in order to take appropriate counter measures. 

The problem is however extremely difficult, because the definition of what is a bubble is 
prone to controversies. Superficially, financial bubbles are easily defined as transient upward 
accelerations of the observed price 

the determination of a bubble requires, in this definition, a precise determination of what 
is the fundamental value. But, the fundamental value is in general poorly constrained. In 
addition, a transient exponential acceleration of the observed price that would be taken as 
the diagnostic of a developing bubble is not distinguishable from an exponentially growing 
fundamental price. 

The Log Periodic Power Law (LPPL) proposed by A. Johansen, O. Ledoit and D. Sornette 
(JLS) 3), 5-7] proposes a way out of this dead-end by defining a bubble as a transient "faster- 
than-exponential" growth, resulting from positive feedbacks. The JLS model provides a 
flexible framework to detect bubbles and predict changes of regime from the study of the 
price time series of a financial asset. In contrast to the more traditional view, the bubble 
is here defined as a faster-than-exponential increase in asset prices, that reflects positive 
feedback loop of higher return anticipations competing with negative feedback spirals of 
crash expectations. The Johansen-Ledoit-Sornette models a bubble price as a power law 
with a finite-time singularity decorated by oscillations with a frequency increasing with 
time. 

The LPPL model in its original form presents a function of 3 linear and 4 nonlinear 
parameters that should be estimated by fitting this function to the log-price time series. 
Calibrating the LPPL model has always been prone with difficulties due to the relatively 
large number of parameters that must be estimated and the strong nonlinear structure of 
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the equation. Most of the fitting procedures that have been used until now subordinate 

ft 

the 3 linear parameters to the 4 nonlinear parameters [6j. The resulting search space with 
4 nonlinear parameters has a very complex quasi-periodic structure with multiple minima. 
The determination of the global minimum requires using some metaheuristic methods such 

nn n 

as taboo search p,|8| or genetic algorithm [9[. But even these methods do not ensure that the 
correct solution is discovered. In addition, how to deal with the existence of many possible 
competing degenerate solutions has not been satisfactorily solved. 

In the present paper, we propose a fundamental revision of the formulation of the LPPL 
model, that transforms it from a function of 3 linear and 4 nonlinear parameters into a 
representation with 4 linear and 3 nonlinear parameters. This transformation significantly 
decreases the complexity of the fitting procedure and improves its stability tremendously 
because the modified cost function is now characterized by good smooth properties with 
in general a single minimum in the case where the model is appropriate to the empirical 
data. As an additional step, we propose a subordination procedure that slaves two of the 
nonlinear parameters to what can be considered to be the most crucial nonlinear parameter, 
the critical time t c defined as the end of the bubble. The critical time is indeed the prize of 
the whole forecasting exercise: the sooner a bubble is identified, that is, the further away 
is t c from the present time, the better it is for policy makers to take appropriate actions. 
Of course, it goes without saying that the validity and reliability of the results should be 
established carefully before any action is taken. We note that this additional subordination 
decreases further the complexity of the search and provides an intuitive representation of the 
results of the calibration. With our proposed methodology, metaheuristics are not longer 
necessary and one can resort solely to rigorous controlled local search algorithms, leading to 
dramatic increase in efficiency 

The paper is organized as follows. Section 2 describe the idea behind the LPPL model 
and presents its original form. Traditional fitting procedures are described in section 3. 
Section 4 presents the proposed modification of the model and explain the new subordination 
procedure. Section 5 concludes. 
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II. LOG-PERIODIC POWER LAW MODEL 



The JLS model [5|— LZJ] assumes that the logarithm of the asset price pit) follows a standard 
diffusive dynamics with varying drift fi(t) in the presence of discrete discontinuous jumps: 

y = /i(t)dt + a(t)dW - ndj . (1) 

In this expression, a(t) is the volatility, dW is the infinitesimal increment of a standard 
Wiener process and dj represents a discontinuous jump such as j = before the crash and 
j = 1 after the crash occurs (// b af f tcr dj = 1). The parameter n quantifies the amplitude of the 
crash when it occurs. The expected value of dj is nothing but the crash hazard rate h(t) times 
in the infinitesimal time increment dt: E[dj] = h(t)dt. The JLS model assumes that two 
types of agents are present in the market: a group of traders with rational expectations and 
a group of noise traders who exhibit herding behavior that may destabilize the asset price. 
According to the JLS model, the actions of noise traders are quantified by the following 
dynamics of the hazard rate [5K7|: 

h(t) = a(t c - t)™- 1 (1 + (3 cos(w ln(t c - t) - </>')), (2) 

where a, (3, u and (j) are parameters and t c is the critical time that corresponds to the end of 
the bubble. The power law behavior (t c — t) m_1 embodies the mechanisms of positive feedback 
at the origin of the formation of bubble. The log-periodic function cos(u ln(t c — t)— <p') takes 
into account the existence of a possible hierarchical cascade of panic acceleration punctuating 
the course of the bubble. The no-arbitrage condition E[dp] = imposes that the excess 
return fi(t) is proportional to the crash hazard rate h(t): fi(t) = nh{t). Solving equation (CQ) 
with fi(t) = Kh(t) and under the condition that no crash has yet occurred (dj = 0) leads to 
the following log-periodic power law (LPPL) equation for the expected value of a log-price: 

E[lnp(t)] = A + B(t c - t) m + C(t c - t) m cos(w ln(t c -t)-<f>), (3) 



where B = —na/m and C = —KaP/^/m 2 + u 2 . It should be noted that solution ([3]) 
describes the dynamics of the average log-price only up to critical time t c and cannot be 
used beyond it. This critical time t c corresponds to the termination of the bubble and 
indicates the change to another regime, which could be a large crash or a change of the 
average growth rate. 
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The LPPL model ([3]) is described by 3 linear parameters (A, B, C) and 4 nonlinear pa- 
rameters (m, u, t c , 0). These parameters are subjected to the following constrains. Since the 
integral of the hazard rate ([2]) over time up to t = t c gives the probability of the occurrence 
of a crash, it should be bounded by 1, which yields the condition m < 1. At the same time, 
the log-price ([3]) should also remain finite for any t < t c , which imply the other condition 
m > 0. In addition, the requirement of the existence of an acceleration of the hazard rate as 
time converges towards t c implies B < 0. Additional constraints emerge from a compilation 
of a significant number of historical bubbles js), Q, 11] that can be summarized as follows: 



0.1<m<0.9, 6<w<13, \C\ < 1, B < 0. (4) 

These conditions (jl]) can be regarded as the "stylized features of LPPL". The condition 
6 < uj < 13 constrains the log-periodic oscillations to be neither too fast (otherwise they 
would fit the random component of the data), nor too slow (otherwise they would provide 



a contribution to the trend (see Ref. 



12| in this respect for the conditions on the statistical 



significance of log-periodicity). The last restriction |C| < 1 in (j3J) was introduced in 13( to 
ensure that the hazard rate h(t) remains always positive. 



III. SUMMARY OF EXISTING FITTING PROCEDURE OF THE LPPL FUNC- 
TION © 

Forecasting the termination of a bubble amounts to finding the best estimation of the 
critical time t c . This requires calibrating the LPPL formula (|3]) to determine the parameters 
t c , m, u, 0, A, B, C of the model that best fit some observed price time series p(t) within 
a time window t G [ti,^]- This may be performed using the Least-Squares Method of 
minimizing the sum of squared residuals: 

N 2 

S(t c , m, u, 0, A, B,C) = J2 [lnp(rO -A- B(t c - n) m - C{t c - n ) m cos(w ln(t c - n) - 0)] , 

i=i 

(5) 

where T\ = t± and — £2- Minimization of such nonlinear multivariate cost function S is 
a non-trivial task due to presence of multiple local minima, where the local optimization 
algorithm can get trapped. However, the complexity of the optimization problem may be 
significantly decreased by noticing that three linear parameters A, B, C can be slaved to the 
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four other nonlinear parameters t c ,m,u,(j) [6j. Indeed, it is easy to prove that 



min S(t c , m, to, 0, A, B, C) = min Si(t c ,m,u, 

' ' ' ' ' t c ,m,uj,ih 



where 



J^hn, S(t c , m, u, (ft, A, B, C). 



Si(t c ,rn,u},(/)] 

The last optimization problem ([7|) may be rewritten as: 



(6) 



(7) 



N 



{A, B, C} = arg min^S(t c , m, u, <p, A, B, C) = arg min | V% — A — Bfa — Cg. 



a,b,c ■ 



(8) 



i=i 



where yi = In pfa), = (t c - Tj) m , g { = (t c - Ti) m cos(u \n(t c - n) - <p). Being linear in 
terms of the variables A, B, C, for m ^ 0, u ^ 0, t c > t 2 and for any (p, this problem has one 
unique solution that an explicit analytical solution obtained from the first order condition. 
The first order condition leads to the matrix equation 



N £/< £<7i 



E/i E/i 2 J2fi9i 
J2fi9i E^ 2 J 



B 



( 



\ 



\ E Vidi J 



(9) 



which is solved in a standard way using the LU decomposition algorithm 14|. Then, the 
global fitting procedure is reduced to solving the nonlinear optimization problem 



{t c , m, u), (J)} = arg min Si(t c ,m,u, 

t c ,m,u},<t> 



(10) 



where the cost function is given by expression (171) . 

Reducing the number of parameters from 7 to 4 simplifies considerably the calibration 
problem. However, the minimization problem f llOp still requires finding the global minimum 
in a four dimensional space of a cost function (J7J) with multiple extrema. The complexity 
of a typical cost function is illustrated by the different cross-sections shown in fig. [2j This 
figure is obtained for the daily time-series of the Shanghai Composite index (SSE) from 
January 2007 to March 2008, shown in fig. HJ Here, for fitting purposes, the time window 
between t x = '12-Mar-2007' and t 2 = '10-Oct-2007' was considered. 

The complex multiple extrema structure of the cost function does not allow one to use 
local search algorithms, such as the steepest descent or the Newton's method, in order to 
find the solution of ([10]) . Such complex optimization problems were studied extensively in 
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FIG. 1: Time dependence of the SSE Composite Index (Shanghai Composite index) from January 
2007 to January 2008 (thin noisy solid line) and the best fit obtained by the LPPL model ([3]) (thick 
solid line). Vertical dashed lines delineate the time window [ti,^] used in the fitting procedure. 

the last 30 years and a number of approaches developed mostly within the class of so-called 
metaheuristic algorithms 15] were introduced to solve this class of problems. Metaheuristic 
algorithms usually make only a few or no assumptions on the cost function, do not require 
calculating the derivatives and can search a very large multidimensional spaces of candidate 
solutions. The cost of such universality and generality is the absence of any guarantee 
of finding the optimum or even a satisfactory near-optimal solution. This issue is usually 
solved with the combination of an initial metaheuristic explorative search followed by a local 
descent to the closest minimum in a second step. 

For the problem of calibrating the LPPL function to financial price time series, the taboo 
search 1(| has been the main approach [g, 0, \v^. Taboo search enhances the performance 
of a local search method by using memory structures that describe the visited solutions: 
once a potential solution has been determined, it is marked as "taboo" area, so that the 
algorithm does not visit explored regions repeatedly. Being metaheuristic, the taboo search 
does not guarantee convergence. Thus in the original work jsj, it was proposed to keep the 
10 best outcomes of the taboo search as initial conditions for a local Levenberg-Marquardt 
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(a)Si(t c ,m= 0.7, u = 7.5, 



(b)5i (t c = '06-Nov-2007', m,u},<f>= tt/10) 




(c)Si(i c ,m = 0.7,u,<t> = tt/10) 



(d)5i(t c = '06-Nov-2007', m, w = 7.5, 



FIG. 2: Different cross-sections of the cost function Si ([7]) for the time-series of SSE Composite 
Index presented at fig. [TJ 



nonlinear least squares algorithm 
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191 ] . The solution with the minimum sum of squares 



between the fitted model and the observations is then taken as the final solution. 

Another important metaheuristic algorithm that was proposed and successfully used to fit 
the LPPL function to financial time series is the genetic algorithm [9]. The genetic algorithm 
mimicks the natural selection process occurring in biological systems, and is governed by 
four phases: a selection mechanism, a breeding mechanism, a mutation mechanism and 
a culling mechanism 20(. Similar to the taboo search, the results of the genetic algorithm 
were refined in a second step by using them as a starting values for the Nelder-Mead simplex 
method 21]. 

Using any of the fitting procedures above with expression always yield some result. 



9 



However, this does not mean that this result should be trusted. The constraints (jl]) con- 
stitute a convenient approach that has proven very useful to filter the solutions that are 
believed to meaningful and relevant to describe developing bubbles from the spurious ones 
(those that correspond to "fitting an elephant" with formulas that contain 4 or more pa- 
rameters. Recall Enrico Fermi as cited by Freeman Dyson who said: "I remember my friend 
Johnny von Neumann used to say, with four parameters I can fit an elephant, and with five 
I can make him wiggle his trunk."). 

As already mentioned, the procedures based on the taboo search and on genetic algo- 
rithms, even when supplemented with the local search algorithms, do not guarantee con- 
vergence to the best solution of the optimization pro blem (I10p . This opens the gates for 

- |25| . In the present paper, we present 



criticisms directed to any of such fitting procedures 22 



a reformulation of the optimization problem f JTQj) that simplifies considerably the fitting 
method. Applied to previous calibrations using the procedures described in this section p- 
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171 . |26[ , we are able to confirm the essential goodness of fits of the obtained calibrations 



on empirical financial time series. 



IV. NEW FITTING METHOD FOR THE LOG-PERIODIC POWER LAW 

This section presents the new calibration method of the LPPL model, which can be 
formulated in two steps that are successively described in the coming two subsections. 



A. From 3 to 4 slaved linear parameters by transformation of the phase 

The key ideas of the new proposed methods is to decrease the number of nonlinear 
parameters and get rid at the same time of the interdependence between the phase <fi and 
the angular log- frequency u. For this, we rewrite the LPPL formula ([3]) by expanding the 
cosine term as follows: 

lnE[p(t)] = A + B(t c -t) m + C(t c -t) m cos(wln(t c -t)) cos 0+C(t c -t) m sm(u\n(t c -t)) sin0. 

(11) 

Let us introduce two new parameters 

Ci = Ccos0, C 2 = Csin0 (12) 
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and rewrite the LPPL equation ([3]) as: 

lnEfo(t)] = A + B(t c -t) m + C 1 (t c -t) m cos{uln{t c -t)) + C 2 (t c ~t) m sm{uln{t c -t)). (13) 

As seen from (fTBl) . the LPPL function has now only 3 nonlinear (t c ,u,m) and 4 linear 
A, B, Ci, C% parameters, and the two new parameters C\ and C 2 contain formerly the phase 



As in the previous section, in order to estimate the parameters, we use the least-squares 
method with cost function 

N 



F(t c ,m,u } A,B,C 1 ,C 2 ) = J2[^p(T i )-A-B(t c -T i y 

i=l 



Ci(t c -r J ) m cos(u;ln(t c -r i ))-C , 2 (t c -r,rsin(u;ln(t c -r i ))J . (14) 

Slaving the 4 linear parameters A, B, Ci, C 2 to the 3 nonlinear t c , u, m, we obtain the non- 
linear optimization problem 



{t c ,m, to} = arg min Fi(t c ,m,uj) 

t c ,m,uj 



where the cost function Fi{t c ,m,uj) is given by 



Fi(t c , m,uj) = ^ jp^ [1 c Fftc, Tn, u, A, B, C%, C2). 



(15) 



(16) 



(17) 



Similarly to the procedure (jHl) leading to (IT71) . the optimization problem ({A, B, Ci, C 2 } = 
argmin J 4 i B i c li c 2 F(t c , m, u>, A, B, C\, C 2 )) has a unique solution obtained from the matrix 
equation: 

/ n E£ T,9i Eh W A \ I J2vi \ 

E/i 12 fi J2fi9i J2fihi B _ J2Vifi 

Y,9i E#A Ci Y.Vi9i 

VE^ E/A E9A E^ 2 / VoJ \J2vihiJ 

where y { = lnp(ii), /; = (t c - T;) m , g t = (t c - T;) m cos(w ln(t c - 7*)) and ^ = (t c - 
Ti) m sin(a; ln(t c - n)). 

The modification from expression to formula (|T3|) leads to two very important results. 

• First, the dimensionality of the nonlinear optimization problem is reduced from a 4- 
dimensional space in (fTOl) to a 3-dimensional space in (jT5l) . This significantly decreases 
the complexity of the problem. 
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• Second, and possibly even more important, the proposed modification eliminates the 
quasi-periodicity of the cost function due to subordination of the phase parameter as 
a part of C\ and C<i to angular log-frequency parameter u. The existence of multiple 
minima of the cost function (as in fig. [2]) has been the main property requiring the use 
of non rigorous metaheuristic searches. The new formulation does not require such 
heuristics and rigorous search methods are now sufficient. 

Fig. [3] presents various cross-sections of the cost function F\ (1161) for the same data set 
that was used in the previous section, namely the SSE Composite Index presented in fig. [TJ 
We use the parameters that were used for fig. [2j One can observe that the cost function now 
enjoys a very smooth structure with only a few (in the fig. [3J — only one) minima, which can 



easily be found using 



squares algorithm 18, 



ocal search methods, such as Levenberg-Marquardt nonlinear least 
19] or the Nelder-Mead simplex method j^ . 



B. Decomposing the optimization search by singularizing the critical time t c 

Being the crucial parameter for forecasting the termination of a bubble, the critical time 
t c requires special care and attention during the fitting procedure. We propose to reformulate 
the optimization problem ( [15]) by using a similar subordination idea that previously allowed 
to separate the linear and nonlinear parameters in (fT5j) , f|T6|) and fflUj) . f|7|) . In this goal and 
without loss of generality, we rewrite problem ( fl5|) as: 

i c = argminF 2 (t c ), (18) 

F 2 {t c ) = minFi(t c ,m,u;), {rh(t c ), u(t c )} = argmin Fi(t c , m, u) (19) 

where Fi(t c ,m,u) is given by ( TT6|) . The optimization procedure (Tl9|) operates on the pa- 
rameter space of the variables m and u. As seen from fig. BK, the cost function Fi(t c , m, u) 
has a very smooth shape and, in the presented case, only one local minimum. Since, in 
the general case, the cost function Fi(t c ,m,u) could have more than one minimum within 
the range of parameters (j3J), we have performed extensive numerical analyses, using the 
SSE Composite index from July 1999 to May 2011 as the empirical case study. We used 
a moving window [ti, ^2] with length of 6 months, scanning the whole range of dates. In 
each window, we counted the number of local minima of Fi(t c , m, u) that fall in the range of 
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(a)F 1 (t c = '06-Nov-2007\m,o;) 



(b)Fi(t CJ m = 0.7 J w) 




(c)Fi(tc,m,a; = 7.5) 

FIG. 3: Different cross-sections of the cost function F\ (|16p for the SSE Composite Index index 
shown in fig. [TJ 

parameters (jlj) for t c varying from t%+l day to £2 + 90 days. The main result is that we never 
found more than three clearly distinguishable local minima of the function Fi(t c ,m,u). In 
other words, the degeneracy in the worst possible cases is very low. We have found that, 
when several minima are present, they can be easily determine by launching no more than 
20 searches with local algorithms started from random points {mo,uo} within the region of 
0.1 < m < 0.9, 6<(jJ <13. We should also stressed that most of the windows that are 
qualified to be in the bubble regime give only one minimum. The occurrence of two to three 
competing minimum has been found to correspond to poor fits of the empirical price time 
series by the LPPL function. 

The cost function i^^c) of the optimization procedure ffl8|) is also very smooth and has 
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only a few minima (see fig. H]for illustration). Again, it is easy to identify these minima by 
using local search algorithms that start from several different initial points t C Q. Such subor- 
dination does not decrease really the computational complexity of the search in fll8p - (TT9|) in 
comparison with the 3-dimensional problem (I15p. In general, the required computations are 
comparable in both approaches. The most important consequence of the reformulation (Il8p 
with ( fl9|) is the possibility of studying the quality of the fit of the critical time t c separately 
and the dependence of the other parameters (m and u) on the critical time: rh(t c ) and (2)(t c ). 
Fig. H] illustrates the dependence of the cost function F 2 {t c ) and of the estimated parameters 
rh(t c ) and oo(t c ) as a function of the critical time t c for the price time series presented in 

fig.m 



T 




Nov07 Dec07 Jan08 Feb08 Mar08 Apr08 



FIG. 4: Plot of the cost function F2(t c ) (top) and of the parameters m(t c ) (middle) and <2>(t c ) 
(bottom) as a function of the critical time t c (defined as the end of the bubble) for the SSE 
Composite Index presented in fig.[TJ The circles correspond to the global minimum {t c , m(t c ),(2)(t c )} 
of the cost function Fi(t c ,m,uj). Dashed lines delineate the domains defined by the constrains ©. 
The thick line represents the values of the parameters where these constrains are met. 
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V. CONCLUDING REMARKS 



This paper provides an important reformulation of the log-periodic power law (LPPL) 
equation of the Johansen-Ledoit-Sornette model of financial bubbles which, in its original 
form, is parametrized by 3 linear and 4 nonlinear parameters. Diagnosing a bubble and 
forecasting its burst (that usually results in dramatic change of regime, such as large crash 
or change of the average growth rate) require calibration of the LPPL model and estimation 
these 7 parameters from the observed price time-series. Though 3 linear parameters could 
be easily slaved to 4 nonlinear ones, the resulting 4-dimensional search space still keeps a 
complex structure with quasi-periodicity and multiple minima. The complex structure of 
the parameter space does not allow one to use local search algorithms and requires more 
sophisticated methods, such as taboo search or genetic algorithms. The implementation and 
the tuning of these methods were the main source of complexity in the calibration of the 
LPPL model. 

The reformulation of the LPPL equation proposed here allows to re-parametrize the 
model in terms of 4 linear and only 3 nonlinear parameters. This has two very important 
consequences, namely: (i) reduce of dimensionality of the nonlinear optimization problem 
from 4-dimensional to 3-dimensional space and (ii) elimination of the quasi-periodicity and 
of the multiple local minima of the cost function. An empirical case study performed in this 
paper using the SSE Composite index from July 1999 to May 2011 shows that, even in the 
worst cases that were studied, we found no more than 3 competing minima. Most cases that 
are qualified to be in a bubble regime give only one minimum, and the occurrence of a second 
or third minimum usually correspond to poor fits of the LPPL function on the empirical 
price time-series. This simplification allows us to use rigorous local search methods in fitting 
procedure without introducing any metaheuristics. 

As an additional step, we have proposed a subordination procedure that slaves two of 
the nonlinear parameters, namely the growth rate exponent m and the log-frequency uj, 
to the critical time t c , that defines the end of the bubble and is considered as the most 
crucial parameter. This additional subordination provides the following benefits: (i) further 
decrease of the complexity and (ii) an intuitive representation, that includes as a bonus an 
explicit dependence of the parameters m and u on the critical time t c . 

Summarizing, the proposed reformulation of the LPPL model decreased dramatically the 
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complexity of the problem and therefore the effort required for its implementation and use 
for the diagnosis and forecasting of bubbles in financial markets. 
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